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Abstract. Dynamic perturbation equations are derived for a generic stationary state 
of an elastic string model - of the kind appropriate for representing a superconduct- 
ing cosmic string - in a flat background. In the case of a circular equilibrium (i.e. 
vorton) state of a closed string loop it is shown that the fundamental axisymmetric 
(n = 0) and lowest order (n = 1) nonaxisymmetric perturbation modes can never be 
unstable. However, stability for modes of higher order (n > 2) is found to be non- 
trivially dependent on the values of the characteristic propagation velocity, c say, of 
longitudinal perturbations and of the corresponding extrinsic perturbation velocity, 
v say. For each mode number the criterion for instability is the existence of nonreal 
roots for a certain cubic eigenvalue equation for the corresponding mode frequency. 
A very simple sufficient but not necessary condition for reality of the roots and there- 
fore absence of instability is that the characteristic velocity ratio, c/vbe greater than 
or equal to unity. Closer examination of the low velocity (experimentally accessi- 
ble) nonrelativistic regime shows that in that limit the criterion for instability is just 
that the dimensionless characteristic ratio c/v be less than a critical value Xc whose 
numerical value is approximately ~- m the relativistic regime that is relevant to su- 
perconducting cosmic strings the situation is rather delicate, calling for more detailed 
investigation that is postponed for future work. 



1. The fundamental dynamical equations of an elastic string model 

The ultimate motivation of the present work is to deal with a problem that was 
implicitly raised when it was first pointed out by Davis and Shellard [1-3] that su- 
perconducting string loops might be "stabilised" (as "vortons") by centrifugal force 
(which can be expected to be very much stronger than the inefficient magnetic 
"spring" repulsion effect that had been previously considered by some authors [4- 
7]) with the implication that (unlike nonconducting cosmic string loops) such current 
carrying loops might survive long enough to produce a cosmological mass excess. 
However, although Davis and Shellard provided the first examples of a general cate- 
gory of circular rotating equilibrium states, their use of the term "stabilisation" was 
rather premature at that stage because they had not carried out any verification of 
the actual stability of such equilibria with respect to ordinary macroscopic dynamical 
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perturbations. (What they did do was to consider the range of parameters for which 
the equilibrium would be able to survive against macroscopic quantum tunnelling.) 
It was in fact fairly easy to confirm that the equilibrium states in question would be 
stable against strictly axisymmetric (circularity conserving) perturbations [8, 9] but 
the present works constitutes what, as far as we know, is the first investigation of 
stability of centrifugally supported string loops against general nonaxisymmetric per- 
turbations. The range of the present study is not limited to the strongly relativistic 
regime relevant to applications to cosmic strings (which will be considerered in more 
detail in a following article) but covers quite general elastic string models including 
those relevant to the experimentally accessible nonrelativistic limit. 

When we first set out on this investigation we had expected that its outcome 
would be expressible simply as a general stability theorem - confirming for generic 
dynamical perturbations what we already knew to be the case for axisymmetric ones. 
However, the picture that actually emerges is less simple and more interesting. It 
turns out that, whereas for many kinds of string model the circular equilibrium states 
will indeed always be stable, nevertheless there are other models for which some or 
even all of the circular equilibrium states are actually unstable. The domain of insta- 
bility is not restricted to the strongly relativistic regime relevant for superconducting 
cosmic string theory and its cosmological implications (which, as remarked above, 
will be considered in more detail in a subsequent article) but it extends down into 
the Newtonian regime. What will be confirmed, however, is that the requirements for 
dynamic stability will always be satisfied in the Hookean (low tension) limit. 

The plan of this work is as follows. After completing this introductory section 
by summarising the general principles governing the dynamics of the category of 
elastic string models under consideration, we go on to provide, in Section 2, a brief 
derivation of the conditions for stationary equilibrium with repect to a nonrotating 
frame in a flat background. Section 3 contains a derivation of the local dynamic 
equations governing the evolution of small perturbations of such equilibrium states, 
ans Section 4 describes the way to obtain global solutions of these equations for cases 
in which the unperturbed equilibrium states are circular (the study of more general 
equilibrium states being left for subsequent work). 

Any two-dimensional worldsheet can be conveniently described [8,10] in terms 
of some chosen orthonormal frame consisting of mutually orthogonal timelike and 
spacelike tangent vectors and as characterised by 

u» Ufl = -l, u% = 0, v% = l, (1.1) 

on the understanding that we are using units such that the speed of light is set to 
unity. In particular, any such frame can be used to provide an explicit representation 
of the antisymmetric tangent element tensor of the surface, which will be expressible 

as 

E^" = vPv v -u v v», (1.2) 

and which is well defined - independently of the choice of frame - modulo an ambiguity 
of sign which is resolved by a choice of orientation in the surface. The square of this 
surface element tensor gives the first fundamental tensor of the worldsheet as 

rf v = gwgv = — u^u v + v»v u , (1.3) 
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which is geometrically well defined, without any ambiguity at all, as also is the second 
fundamental tensor K f, which is defined (for an embedded surface of arbitrary 
dimension) in terms of the corresponding first fundamental tensor by the formula 

K ll S = r h l r VvV* p , (1-4) 
where the surface tangential covariant differentiation operator is defined by 

% = vSVv (1-5) 
As well as the obvious orthogonality and tangentiality properties 

iW = 0, _L^K M / = 0, (1.6) 

where the orthogonal projection operator is defined by 

-11= (1-7) 

the second fundamental tensor has the nontrivial Weingarten symmetry property [8, 
11, 12] expressible - using square brackets to denote index antisymmetrisation - as 

K [liV f = 0. (1.8) 

As in the preceding investigation of local stability conditions and characteristic 
speeds [8, 11], condition (1.8) will play an important role in the calculation that 
follows, since although it will hold automatically as an identity if K^f is derived 
from a two-surface that is already specified, on the other hand it must be deliberately 
imposed as an indispensible (Frobenius type) integrability condition if one chooses, 
as we shall do in this work, to use an approach in which the surface is considered as a 
secondary construct determined implicitly by prior specification of the surface element 
E^ v . Our reason for using such an approach - i.e. for considering the world sheet 
to be constructed from the tangent tensor £^ v rather than the other way round - is 
that, although our ultimate purpos is to consider the world sheet of a single isolated 
string loop, it will be computationally convenient to consider it as a member of a 
fictitious space-filling congruence of string loops, each separately subject to the same 
dynamical laws. This will enable us to treat quantities such as £^ u as ordinary space- 
filling tensor fields free from the restriction of having their support restricted to any 
single world sheet. This facilitates the use of an economically covariant formalism 
that avoids the need to become involved in the questions of choice of gauge that 
would arise if we used an approach of the more traditional kind involving the explicit 
specification of the string worldsheet in some particular coordinate gauge. 

In addition to the purely geometrical evolution condition (1.8) the calculation 
that follows will depend on the use of appropriate extrinsic and internal dynamical 
equations. For any classical string model - in the limit in which finite thickness effects 
including polarisation are neglected - the extrinsic dynamical equation of motion 
governing the motion of the world sheet will be given, whenever external forces are 
absent as is postulated in the present work, by a "generalised sail equation" [8, 11, 
12] of the form 

T^K M / = 0, (1.9) 
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where T^ v is the relevant (surface) stress momentum energy density tensor, whose 
specification will depend on the kind of string model under consideration. (For ex- 
ample, in the well known degenerate cas of a Goto-Nambu string model the tensor 
T^ u is simply proportional to the fundamental tensor r\^ v . 

There is no loss of generality in considering T^ v to be specified by its eigenvalues, 
T and U say, which are interpretable respectively as the intrinsic string tension and 
as the energy density in the preferred rest frame determined by its spacelike and 
timelike eigenvalues v^ 1 and w M which can conveniently be used to specify the choice 
of the frame introduced in (1.1). The surface stress momentum energy density is thus 
expressible in the generic form 

fr 1 " = Uu"u v -Tv^v" , (1.10) 

where what distinguishes one kind of string model from another are the rules speci- 
fying the energy density U and the tension T. 

Although more complicated possibilities can be envisaged (and would be neces- 
sary for dealing with cases in which internal dissipation is important), our present 
study will be limited to the category of simple elastic string models, in which the 
rules for specifying U and T consist just of a single equation of state specifying one of 
them as a monotonic function of the other. (The Goto-Nambu case can be classified 
within this category as the degenerate limit case in which U and T are fixed with 
the same value U = T = m 2 , where m is a constant with the dimensionality of mass 
in units of the kind used in the present work for which the speed of light is set to 
unity.) As well as providing a good description for such familiar examples as ordi- 
nary violin strings, such simple elastic string models are also [8, 9, 13] appropriate 
for the more exotic case of superconducting cosmic string models whose investigation 
provided the primary motivation for this work and which are adequately described by 
the formalism used here, provided that electromagnetic coupling corrections (of the 
kind investigated by Peter [14, 15] are unimportant - as will commonly be the case 
in the cosmological applications we have in mind [1-3, 16, 17]. 

Whatever the equation of state may be, the internal equations of motion for such 
a simple elastic string model can be shown [8, 9] to reduce just to a pair of surface 
current conservation laws expressible as 

V>O=0, (1.11) 

and 

V P OuO=0, (1.12) 

where v and fi are a number density variable and an associated effective mass vari- 
able - or chemical potential - which are determined by the equation of state. It is 
convenient to consider the equation of state to be specified by the expression for the 
energy density U as a function of the number density v. The corresponding function 
for the chemical potential fi will then be given as the derivative 

fi = dU/du, (1.13) 

while finally the function for the string tension T will be given by 

T = U-Pfi. (1.14] 
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It is to be remarked for reference later on that the number consrvation law (1.11) 
can be expressed in the equivalent adjoint form 

£ pa V p (uv a ) = 0, (1.15) 

while the analogous irrotationality condition obtained as the adjoint of (1.12) is 

£ pt7 V p (fiu a ) =0. (1.15) 

In terms of the preferred reference frame introduced by (1.10) the extrinsic equation 
of motion (1.9) for the world sheet can be written in the more explicit form 

±% {Uu v V u u p - Tv"V v v p ) = , (1.17) 

while, finally, the geometric integrability condition (1.8) will take the form 

_L/* p (u v V v v p - v v V v u p ) = . (1.18) 



2. Equilibrium with respect to a static Killing vector in a flat background 

The concept of equilibrium is generally understood to mean invariance with re- 
spect to an at least locally timelike Killing vector /c M of the relevant background 
spacetime, where the Killing condition for the invariance of the background space- 
time metric g^ v itself is expressible in terms of Riemannian covariant differentiation 

as 

V (/ A)=0, (2.1) 

using round brackets to indicate symmetrisation. The Killing vector is said to be 
static in the special case in which it also satisfies the irrotationality condition 

V^A>]=°- ( 2 - 2 ) 

In the case of a flat background spacetime, the combination of the Killing condition 
(2.1) and the staticity condition (2.2) is expressible simply as the condition of being 
covariantly constant, i.e. 

V M /^ = 0, (2.3) 

which is the equation characterising a simple linear translation. Subject to the appro- 
priate signature convention, we assure that the generator k^ is timelike by imposing 
the usual normalisation condition 

fc% = -1 . (2.4) 

The requirement that the string should be in equilkibrium with respect to k^ 
includes as a first - purely topological - condition the requirement that /c M should be 
tangential to its world sheet. It follows that at each point in the world sheet it will be 
possible to choose an orthogonal spacelike unit tangent vector say, as characterised 
by 

e% = l, e*% = 0, (2.5) 
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in terms of which the fundamental tensor of the world sheet will be expressible as 

rf v = -k»k v + e»e v , (2.6) 

where e M is uniquely determined by a single global choice of orientation on the (topo- 
logically cylindrical) world sheet of the string loop. 

The condiion of equilibrium with respect to a vector field k p means that any 
physically well defined tensor field on the world sheet should be invariant under its 
action. With respect to a vector k p satisfying the covariant constancy condition (2.3) 
(which implies that the corresponding Lie differentation operation reduces to one of 
direct covariant differentiation) the condition of invariance simply means that the 
corresponding covariant derivative should vanish. Thus in the case of the spacelike 
unit vector e M the invariance condition is simply 

^" = 0. (2.7) 

Using this in conjunction with (2.3) we can immediately evaluate the second funda- 
mental tensor (1.4) of the world sheet as 

K^f = e^K? , (2.8) 

where its trace, the curvature vector K p , is given by 

K p = Ky = e u V„e p . (2.9) 

Although we are restricting ourselves to cases of equilibrium with respect to a 
Killing vector that is static (i.e. the present work does not include consideration 
of equilibrium with respect to a rotating frame) we cannot impose the requirement 
that the equilibrium itself be static in the strictest sense - which would mean that 
the intrinsically preferred unit vector u p of the string would be aligned with the 
Killing vector generator k^ - since (as shown by previous more general studies [18, 
19]) such a strictly static equilibrium is possible only for states in which the string 
is actually straight, meaning that K^ 1 should vanish everywhere along it, a condition 
which is obviously incompatible with the closed loop topology under consideration 
here. The kind of equilibrium states to be investigated here are not strictly static 
but merely stationary, which means they will be characterised by a nonzero "running 
velocity", v say, expressing the deviation of the intrinsically preferred timelike unit 
stress momentum energy eigenvector of the string from the Killing vector so 
that it will be expressible in the form 

u p = 1 (k p + ve fI ) , (2.10) 

while similarly (subject to the appropriate sign convention) the corresponding space- 
like unit stress momentum energy eigenvector v M will be expressible in the form 

v^ = 1 (vk^ + e^), (2.11) 
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where 7 is the relevant Lorentz factor, i.e. 



7 = 1/^1-^2. (2.12) 

It follows from the results of the previous paragraph that the derivatives appearing 
in the integrability condition (1.18) will be expressible as 

u u V v v^ = v v V v u» = v<y 2 K^ , (2.13) 

while those appearing in the equation of extrinsic motion (1.17) will be given by 

u v V v v? = v 2 ~f 2 K^ , v u V„v» = ~f 2 K^ . (2.14) 

Substituting this in (1.17) - or by direct substitution of (2.8) in (1.9) - it can be seen 
that the extrinsic dynamical equation will reduce simply to 

(Uv 2 -T)K» = 0. (2.15) 

In so far as the longitudinal dynamical equations are concerned, it can be seen that the 
combination of the irrotationality condition (1.16) with the equilibrium requirement 
that derivatives with repect to should vanish leads to a conservation law of the 
form 

V p ( 7 ^)=0, (2.16) 

which is the string analogue of the well known Bernoulli conservation law in irro- 
tational perfect fluid theory. In the string case (unlike the perfect fluid case) one 
can also derive a corresponding dual conservation law from the other longitudinal 
dynamical equation (1.15), which similarly leads to 

V p (v>yu) = 0. (2.17) 

For a generic equation of state (the noteworthy exception being the special case of the 
integrable equation of state characterised by the condition that the product of U and 
T be constant [8, 19, 21]) the constancy conditions (2.15) and (2.17) are independent 
and therefore sufficient to determine the constancy of the intrinsic state of the string, 
as specified by v, i.e. 

V p z/ = 0. (2.18) 

and also of the running velocity v. 

It can be seen from the extrinsic equilibrium condition (2.15) that there are two 
qualitatively different ways in which local equilibrium with respect to a static Killing 
vector in a flat background can be achieved. One possibiity is for the string to be 
straight, = 0, in which case the running velocity v can have any (subliminal) 
value, including the zero value that correponds to the strictly static case. The other 
possibility - the only one that is relevant when the curvature K p is nonzero as must 
be the case for the closed loops with which the present work is concerned - is that of 
stationary but nonstatic states in which the running velocity has the nonzero value 
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determined uniquely (modulo a sign) as a function of the intrinsic string state by the 
formula 

v 2 = T/U, (2.19) 

(which must hold not only in the generic case for which, according to (2.18) both 
the running velocity v and the intrinsic string state, as characterised by T and U are 
constant, but also in the exceptional case of the integrable string model[19] for which 
the running velocity v and the intrinsic state may vary along the length of the string). 

The formula (2.19) is the same as the formula [8, 11] for the characteristic speed 
(relative to the intrinsically preferred rest frame specified by w M ) of extrinsic perturba- 
tions of the world sheet. The result that has just been derived is therefore expressible 
as the theorem [18, 19] that in any curved equilibrium state of a simple elastic string 
model with respect to a static (non-rotating) Killing vector in a flat background the 
running velocity must be just such as to cancel the characteristic speed of relatively 
backward moving extrinsic perturbations, which are therefore static in the sense that 
the relative velocity, V- say, of such relatively backward moving perturbations with 
respect to the static Killing vector /c M is zero, i.e. 

v_=0. 2.20 

It follows that, with respect to the speed, v + say, of relatively forward moving 
extrinsic perturbations will be approximately twice the relative runing speed v in the 
Newtonian limit in which the relative speeds are small compared with the speed of 
light (which is unity in the system we are using), the exact expression that will be 
valid for arbitrarily high velocities being given by 

2v 

v+ = TT5 5 ' (2 ' 21) 

We remark that the relative characteristic speed, c say, of longitudinal perturba- 
tions, which is given [8, 11] by 

o dT v du 

may, depending on the equation of state, be less than or greater than or equal to 
that of the extrinsic (transverse) perturbations as given by the formula (2.19) for v 
and that in the former case the relatively backward moving longitudinal perturbations 
would be absolutely forward moving relative to the background frame specified by k^. 
Both from the point of view of secular stability with respect to radiation backreaction 
(whose analysis, to which we have alluded elsewhere [12, 19] is beyond the scope of 
the present discussion) and, as will be seen below, from the point of view of ordinary 
dynamical stability, it is the latter alternative c > v, that is most favorable. As 
well as applying to the familiar nonrelativistic Hookean (low tension) limit, the strict 
inequality c > v will hold in particular for an equation of state of the simple constant 
trace type (i.e. a constant value for the sum U + T) that was implicitly used in 
many of the early discussions of superconducting cosmic strings [22-25], for which 
the longitudinal characteristic speed is actually equal to that of light, i.e. c = 1. 
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However, closer examination of the underlying superconducting string teory by Peter 
[13] has shown that it would be more accurate to use an equation of state for which 
the longitudinal characteristic speed is in the dangerous regime, meaning that it is 
lower than the transverse characteristic speed, i.e. c < v, while the special integrable 
string model referred to above [20, 21] corresponds to te intermediate case c = v. 

3. Local perturbation equations for acurved equilibrium state. 

Before proceeding, we need to introduce a mathematical device that is frequently 
useful [11] for dealing with dimensionally restricted physical models - such as the 
string case under consideration here - when one wishes to avoid getting involved in 
the technicalities of Lagragian perturbation theory. In order to be able to work in 
terms of ordinary Eulerian perturbations, i.e. direct comparisons between perturbed 
and unperturbed configurations at the same fixed point in the background spacetime, 
one needs to use a model that spreads out beyond the single "physical" worldsheet 
in which one is ultimately interested. What is required for our purpose is a locally 
extended mathematical model constructed in terms not just of the single "physical" 
world sheet with which one is ultimately concerned, but also of a locally surrounding 
spacefilling congruence of "fictitious" world sheets, each of which is postulated to be 
separately subject to the same mathematical equations that are obeyed by the central 
"physical" world sheet within it. 

Since it can be seen from (2.15) that the space configuration of an equilibrium con- 
figuration is arbitrary provided (2.19) is satisfied, it is easy to extend any "physical" 
equilibrium configuation of a string loop to a neighbouring congruence of equilibrium 
configurations simply by extending the field by parallel propagation orthogonal 
to the central "physical' world sheet, which can be seen from (2.9) to mean that its 
covariant derivative at any point on the physical worldsheet will be given in terms of 
the curvature vector of the latter by 

V li e v = e li K v . (3.1) 

Since this construction repects the stationarity condition (2.7) the necessary integra- 
bility conditions are satisfied as a triviality: the field e M set up in this way automati- 
cally commutes with the Killing vector field so that the two vector fields together 
will necessarily generate a well defined congruence of stationary world sheets. To 
complete the construction of the "unperturbed" state of the fictitious string congru- 
ence we postulate that the equilibrium conditions (2.16) to (2.19) are satisfied on each 
separate world sheet. 

Under such conditions all that is needed for the equilibrium state of the extended 
model to be specified completely is the choice of the running velocity v on each sepa- 
rate world sheet. The simplest possibility, consistent with the logitudinal uniformity 
condition (2.18) is to take the running velocity and, hence, the intrinsic state of each 
string (as specified by its tension T, for example) to be uniform throughout. How- 
ever, in order to retain the freedom of making alternative choices (for the purpose of 
checking consistency as a precaution against algebraic error) we shall provisionally 
retain the option of letting v vary from one world sheet to another as an arbitrary 
differentiable field. 
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Having chosen an equilibrium state in the manner just described, we now consider 
the dynamical evolution of an arbitrary small perturbation of it. Using the symbol 5 
to indicate an Eulerian differential, we can write the perturbations of the longitudinal 
dynamical equations (1.15) and (1.16) directly as 

S(S^)V p (uv a ) + £ pCT V p 8{vv a ) = , (3.2) 

and 

V P 0m ff ) + E^VpSi/jLv,,) = . (3.3) 

After simplification by the use of the unperturbed equilibrium equations the corre- 
sponding perturbation of the trensverse dynamical equation (1.17) takes the form 

{8(Uu v )W u v p - 8{Tv v )V v u p + Uu v V„8{v p ) - Tv v VJ{u p ) . (3.4) 

Finally as the (very modest) price to be paid for the convenience of working in terms 
of simple Eulerian perturbations of unrestricted fields, we must complete the dynamic 
system by including the geometric surface integrability condition (1.8), whose satisfac- 
tion in this approach is not guaranteed automatically as it would be in a Lagrangian 
treatment. Taking the direct perturbation of the explicit version (1.18) of (1.8), and 
simplifying by the use of the equilibrium equations in the same way as was done for 
(3.4), one obtains the corresponding geometric surface integrability condition for the 
perturbed field variables in the form 

{8{u v )V u v p - 8{v v )V v u p + u u W v 8(v p ) - v v V v 8{u p ) . (3.5) 

In order to proceed it is convenient to represent the longitudinal part of the 
perturbation in terms of a pair of infinitesimal scalars, a and e say, and to represent 
the transverse part of the perturbation in terms of a pair of vectors, and say, 
subject to the orthogonality conditions 

rf v & = , (3.6) 

according to the following further specifications. Subject to the normalisation condi- 
tions (1.1) the most general perturbation of the frame vectors -u M and is determined 
by and determines the vectors £± and the scalar e by the relations 

Su" =ev» + v(£ +C+), Sv" =eu»+v(C ~ (+) , (3.7) 

where the coefficient v is the extrinsic propagation velocity as given by (2.19) which is 
included in order for and to be interpretable as the respectively "left" (relatively 
backward) and "right" (relatively forward) moving parts of the extrinsic perturbation. 
Finally the other infinitesimal amplitude a is defined as the logarithmic variation of 
the number density v that was chosen in Section 1 as the basic independent variable 
for characterising the intrinsic equation of state, i.e. we set 

^ = a. (3.8) 
v 
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The corresponding formula for the variation of the chemical potential or effective mass 
variable fx that is required for the evaluation of (3.3) will then be obtainable using 
the equation of state as 

— = c a , (3.9) 

where c is the longitudinal characteristic speed as given by the formula (2.22). As 
remarked above, in the constant trace model used in most early and many more recent 
discussions [22-25] of cosmic string dynamics the longitudinal perturbation speed is 
that of light, i.e. c = 1. In the more accurate models derived by Peter [13] its value is 
not just less than that of light but even less than the transverse characteristic speed, 
i.e. one has c < v < 1. On the other hand, in the low tension limit that has been 
experimentally familiar in the nonrelativistic regime since the time of Newton and 
Hooke, one has v<c<l. 

Using condition (3.1) and the fact that the orthogonality property (3.6) enables 
the contracted surface derivatives of the extrinsic perturbation vectors to be expressed 
in the form 

V M C£ = ~K»e± , (3.10) 

we can rewrite the pair of longituninal perturbation equation (3.2) and (3.4) in char- 
acteristic form, in terms of the quantities that have just been introduced, as the pair 
of equations 

K±c^)V P (£±ca) = ^^{(v±c)C + (vtc)C p + }V p u + 2v 1 2 (v±c)K p C p + . (3.11) 

The pair of extrinsic equations can similarly be combined as a forward characteristic 
equation, 

2v ± a p (u v + vv v )V v C + = -K a {2ve + (v 2 + c 2 )a)} , (3.12) 

in which it may be remarked that does not appear at all, together with the inte- 
grability condition itself which takes the form 

K - w»)V„£ =±% K + iw")V„C£ • (3.13) 

It can be seen that the longitudinal equations (3.11) will greatly simplify if we now 
fix the choice of the unperturbed "fictitious" string states of the extended system by 
requiring that the intrinsic state as characterised by v should be uniform throughout, 
i.e. 

V p z/ = 0, (3.14) 

not just longitudinally along each individual world sheet as required by (2.18). In 
that case (3.11) and (3.12) will form a decoupled system of equations for the two 
longitudinal variables a and e and the two independent components of the foreward 
moving part of the extrinsic perturbation. The remaining equation (3.13) is the 
only one involving the other two independent components, namely those of the rel- 
atively backward mpving part of the extrinsic perturbation. Application of the 
theorem mentionned at the end of Section 2 to the effect that the relatively back- 
ward moving extrinsic perturbation should have zero propagation speed relative to 
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the static Killing vector can be seen directly from the possibility of rewriting (3.13) 
in the form 

-L CT , WVvCt =±% 7 2 (1 + v 2 W + v + enV„( p + , (3.15) 

where v + is the foreward extrinsic propagation speed given by (2.21). The decoupled 
system for the other four independent components is given by the analogous version 
of the characteristic evolution equation (3.12) for the foreward moving part of the 
extrinsic perturbation, i.e. 

2 7 v(l + v 2 ) r p (k v + v + e v )V u C + = -K a {2ve + (v 2 + c 2 )a} , (3.16) 

together with the longitudinal equations which are expressible as 

(1 + v 2 )(k v + v+e v )V„s + {v(l + c 2 )k v + (v 2 + c 2 )e v }V„a = A 1V 2 K P C P + , (3.17) 

and 

k v V v e - 7 2 {f(l - c 2 )k v + (v 2 - c 2 )e v }V v a = . (3.18) 



4. Global solution of the perturbation equations in the circular case. 

In order to obtain a system whose global solution is easily calculable, we now 
restrict ourselves to the case in which the unperturbed equilibrium state is a circular 
ring configuration with fixed radius, r say, which means that the tangential derivative 
of the curvature vector will be given by 

ri%V v K^ = -K v K v e p e^ (4.1) 

where its magnitude has the fixed value given by 

K»K^ = 1/r 2 . (4.2) 

For such a configuration the perturbation may be conveniently analysed in terms 
of the natural time and angle coordinates which are specified (modulo a choice of 
origin) on the equilibrium world sheet by 

TO = 1, W = 0, e u Va = 0, e v V v e = l/r. (4.3) 

We can consider a quite general perturbation as a discrete sum of Fourier components 
with respect to the periodic coordinate 9. For the scalar longitudinal perturbation 
variables a and e the Fourier components are definable for each (by convention pos- 
itive) integer n = 0, 1, 2, simply as the real parts of corresponding complex variables 
having the form 

a = ae- ind , e = ee~ ine , (4.4) 

where a and s satisfy the rotational invariance condition 

e u V u a = 0, e u V u e = 0, (4.5) 
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so that they are functions only of t but not of 9. The corresponding Fourier com- 
ponents of the extrinsic displacement vectors £± will be similarly expressible in the 
form 

C£ = e~™ e , (4.6) 

where the complex amplitude vectors C± are rotationally invariant in the appropriate 
sense, which is not that of covariant differentiation but that of Lie differentiation with 
respect to the relevant angle Killing vector, re M . This Lie type constancy condition is 
just the natural generalisation of the condition (4.1), satisfied by the curvature vector 
itself, being expressible in terms of straightforeward covariant differentiation as 

r? p ^a = -K v C ± e p e». (4.7) 

Condition (2.3) ensures that with respect to the time Killing vector /c M there 
will be no distinction between covariant and Lie differentiation which we shall denote 
simply by a dot. Thus for any Fourier mode of the kind just described, it can be seen 
that the derivations in the perturbation equations obtained in the previous section 
can be evaluated using the simple scalar expressions, 

Tl Tl 

k"V„a = a, k v V u e = e, e u V v a = -i-a, e"V v s = -i-e, (4.8) 

r r 

together with the corresponding vectorial expressions, 

77 

fc"v„C£ = C£, ±» p e»v u ( p ± = -i-(£. (4.9) 

Substitution of these relations enables us to reduce Eq. (3.15) for the relatively 
backward moving part of the extrinsic perturbation to the ordinary differential 
relation 

(i-v 2 )c = (i+v 2 )e + -^e + . (4.io) 

The decoupled system (3.16) - (3.18) for the remaining variables a, e, and the 
two independent components of is analogously reducible to the ordinary differential 

system, 

2 7 v((l + v 2 )(» + ^<#) = (2ve + (v 2 + c 2 )a) , (4.11) 

(1 + J)e - 2 -^l £ + v(l + c 2 )a - (v 2 + c 2 )-a = 2(3 , (4.12) 
r r 

(1 - v 2 )e - v{\ - c 2 )a + (v 2 - c 2 )— a = O , (4.13) 

r 

using the abbreviation 

p = 2 1V 2 K u (+- (4.14) 

The generic solutionof the foregoing evolution equations will have the form of a 
normal mode as given in terms of some constant real or complex u> by 

& = (a)e* at , e=(e)e i " t , g = <C£> e*"* , (4-15) 
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where the (complex) bracket labelled quantities are all constant - not just with respect 
to the angle but also with respect to time, i.e. 



(d) = 0, («) = 0, (C£> = 0. (4.16) 

It can be seen from (4.10) that, as was to be expected from the general equi- 
librium theory described in section 2, for all values of n there are effectively static 
(zero eigenfrequency) perturbation modes for arbitrary time independent complex 
amplitudes of the relatively backward moving extrinsic displacement vector with 
vanishing values of the relatively foreward part and of the longitudinal perturbation 
amplitudes, i.e. modes characterised by 

c^ = 0, C = (C)e- ind C^=0 a = e = 0. (4.17) 

It can also be seen that the azimuthal component of as projected orthogonally 
to the curvature (i.e. parallel to the axis of symmetry of the ring) will decouple 
completely from the rest of the system, which shows that for each value of n there will 
be a corresponding relatively foreward propagating extrinsic mode with displacement 
parallel to the axis of symmetry that takes the normal mode form (4.15) with 

uj = nQ + , C-=0, (3 = a = e = 0, (4.18) 

using a notation scheme in which Q and 0+ denote the angular velocities correspond- 
ing to the unperturbed ring rotation velocity v as given by (2.19) and the relatively 
foreward rotating extrinsic velocity v + as given by (2.21), i.e. 

O = - , 0+ = ^± , 4.19) 

r 1 r 

so that (4.18) can be interpreted as meaning that relatively foreward moving az- 
imuthal modes consisting purely of extrinsic displacements parallel to the axis of 
symmetry will propagate with the same (characteristic) velocity v + independently of 
their wavelength. 

As the only nontrivial (wavelength dependent) part of the proble we are thus 
left with the subsystem of longitudinal and purely radial perturbations as specified in 
terms of just three perturbation variables which can be listed as the scalar (3 defined to 
be proportional to the radial component of the perturbation vector by the relation 
(4.14), together with the longitudinal perturbation amplitudes e and a. The system 
governing the normal modes 

e = ( £ ) e *M-^) , a = (a) e*^*"^) , (3 = ((3) e*^*-**) , (4.20) 

of these three remaining variables is expressible in the matrix form 

2 v 2 + c 2 (l+v 2 )a-' 

(l+v 2 )a-2n v 2 (c 2 + l)a-(v 2 +c 2 )n 2 | ( a | = j j (4.21) 

(l-v 2 )a v 2 (c 2 -l)a + (v 2 -c 2 )n 













a 








[i(3r J 
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in which the perturbation frequency variable u> has been replaced by the corresponding 
dimensionless ratio 

a = u/Vt . 

The system (4.21) evidently provides an eigenvalue equation for the frequency variable 
a, which must be such that the corresponding determinant vanishes, i.e. 



v 2 (l+v 2 )(l-c 2 v 2 )a 3 + 2v 2 (c 2 -v 2 -2(1- c 2 v 2 ))na 2 

+{Av 2 {l -c 2 ) (n 2 - 1) - (1 + v 2 ) (c 2 - v 2 ) (n 2 + 1)) a + 2(c 2 -v 2 ) (n 2 - l)n = . 



(4.23) 



The condition for strict dynamic stability of the mode specified by a given value 
of the integer n is that the corresponding eigenvalue equation (4.23) should have three 
real roots, distinct from each other - and from the zero eigenvalue u> = corresponding 
to the solution (4.17) - in which case there will exist corresponding real values for 
the triplet of amplitudes (e), (a), and i((3). The presence of the factor % in the last 
of these is interpretable as meaning that there is a phase lag whereby the maximum 
variation of the intrinsic variables (such as the tension T) as determined by a, and 
of the longitudinal velocity variation, as determined by e, will always coincide with 
the phase at which the radial displacement as determined by (3 is zero. The existence 
of complex roots will evidently imply strict dynamical instability, i.e. a mode that 
grows exponentially, while the special case of zero or mutually coincident eigenvalue 
solutions of (4.23) corresponds to the existence of a mode that is only marginally 
stable, not strictly stable, so that the normal mode form (4.20) may not be sufficient 
to provide the complete solution of its time evolutionas in the generic case. It is 
evident that (4.23) has a zero eigensolution only for the trivial case n = of a nearby 
concentric circular equilibrium state and for the almost equally trivial case n = 1, for 
which there are marginally stable state that grow not exponentially but linearly in 
time, and that are interpretable merely as equilibrium states with respect to nearby, 
slightly Lorentz boosted frames. 

5. Demonstrations of stability in particular cases. 

It is easy to check that, subject to the standard conditions [8, 11] for microscopic 
stability and causality, i.e. 

0<v 2 <l, 0<c 2 <l, (5.1) 

the conditions for dynamical stability will always be satisfied for the lowest modes 
that are characterised by n = and n = 1. 

In the case n = 0, corresponding to modes of pure conformal expansion or con- 
traction in which the ring conserves its circular shape, the odd powers of a in (4.23) 
drop out, so that the equation is easily solved to give either the trivial possibility 
a = - corresponding to a perturbation in which the original ring solution is merely 
replaced by another neighbouring ring state - or else the fundamental oscilation fre- 
quency which can easily be seen to be given by 



c 2 - 3cV + 3v 2 



v 



4 



v 2 (l + v 2 )(l - C 2 V 2 ) 
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In the first general study of circular equilibrium states of the kind we are considering 
it was argued that stability with respect to axisymmetry preserving perturbatins 
could be derived from the condition that the relevant total mass energy M. should be 
a minimum. That requirement can be seen to be precisely equivalent to the reality 
requirement that must be satisfied by the dimensionless frequency a since the formula 
that was obtained [8, 9] for the second derivative of M. can be seen to be related to 
the fundamental angular frequency u> = fia obtained from (5.2) by the very simple 
relation 

2 1 d 2 M 

U= M^- (5 - 3) 

It was remarked in the original discussion that the right hand side of (5.3) would be 
positive, ensuring stability of the fundamental mode, in all the most obvious appli- 
cations. A closer scrutiny shows that this is actually an understatement: the right 
hand side of (5.2) will in fact be positive in all cases without exception, provided the 
inequalities (5.1) are not violated, as can be seen by regrouping the terms in the form 

2 _ 2v 2 (l-c 2 ) + (c 2 + v 2 )(l-v 2 ) 
a v 2 (l + v 2 )(l-c 2 v 2 ) ' 1 ' 

When we go on to consider the next level of excitation as characterised by n = 1 
we find that since the inhomogenious term in (4.23) still drops out we are again left 
with a trivial zero eigenfrequency (which corresponds to the possibility of a simple 
lateral displacement and whose degeneracy with (4.17) signals the existence of the 
marginally stable linearly growing mode corresponding to the further trivial possibility 
of a simple Lorentz boost), together with the pair of non-trivial oscillation frequencies 
that are obtained as solutions of the quadratic equation 

v 2 (l+v 2 )(l-c 2 v 2 )a 2 + 2v 2 (c 2 -v 2 -2(l-c 2 v 2 ))a - 2(c 2 -v 2 )(l+v 2 ) = . (5.5) 

By casting the explicit solution in the form 

v(v 2 -c 2 + 2(l- c 2 v 2 )) ± ^/c 2 {l-v 4 )+v 2 {l-c 2 v 2 )^{l-v 4 ) + (l-c 2 v 2 ) 
°~ v(l + v 2 )(l - c 2 v 2 ) (5-6) 

it can be seen that, as in the previous case, the roots will always be real so long as 
the inequalities (5.1) are satisfied, since all the bracketted combinations within the 
square roots are manifestly positive. 

For the purpose of investigating higher modes it is useful to use a rescaled variable 

* = i±-%, (5.7) 

and to work in terms of a dimensionless "distension" parameter, 

c 2 - v 2 

A = ^n 2"2V ( 5 - 8 ) 

f+(l — c 2 v 2 ) 

which can be seen to play a critical discriminating role. In terms of this parameter 
A and of the prograde rotation speed v + as given by (2.21) instead of the original 
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parameters c and v, the eigenvalue equation (4.23) can be rewritten in the convenient 
form 

(*-l-^)(*-l^)=A(l-^)(*-l + ; i)^, (5.9) 

whose solution can be interpreted as the geometric problem of finding the intersections 
of the parameter independent parabola determined by the left hand side with the 
hyperbola determined by the righthand side. Noting that the vertical asymptote 
to the hyperbola occurs where x = —A, it can easily be seen that there will be a 
real negative root (corresponding to a retrograde mode) if and only if the distension 
parameter A is positive, and that if this condition is satisfied there will also be two 
real positive roots (corresponding to prograde modes) . The condition of positivity of 
A is evidently equivalent to the condition that the dimensionless ratio 

X = - (5.10) 

v 

of longitudinal to extrinsic characteristic speed should be greater that unity, 

X>1. (5.11) 

This sufficient condition for stability can obviously be extended to cover the case 
for which the discriminant vanishes, i.e. for which 

X = l. 5.12) 

This special condition automatically holds in the special case of the integrable "con- 
stant product" model with equation of state UT = m 4 for some fixed mass parameter 
m - which arises naturally from Kaluza Klein theory [8] and is also relevant as a 
macroscopic (effectively renormalised) "warm" string approximation allowing [20, 21] 
for a thermal or other distribution of microscopic wiggles on an underlying Goto - 
Nambu model - and which is characterised by the condition that the two kinds of 
characteristic speed remain exactly equal as the tension varies. Since this property al- 
lows the equations of motion to be easily integrated explicitly in a flat background [20] 
it could have been predicted in advance that the present perturbation problem should 
also be easily soluble in this particular case. Such an expectation is indeed justified, 
since when the condition (5.12) holds the constant term (4.23) evidently drops out so 
that again, as for the other cases n = and n = 1 one obtains a zero eigenfrequency 
uj = while the other two eigenfrequencies are got by solving a quadratic which leads 
to a result that is expressible very simply as 

u=(n±l)n+, (5.13) 

where 0_|_ is the angular velocity of relatively foreward characteristic motion as given 
by (4.19). The manifest reality of this result confirms that for this particular model 
no instability occurs for any mode. 

Knowing that the conditions (5.1) ensure stability in the large n limit [11] and 
having found no possibility of instability in the lowest modes, n = and n = 1, nor 
for any value of n when % > 1, one might be tempted to conjecture that all modes will 
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necessarily be stable. To show that, despite the existence of all these positive results, 
unstable counterexamples can occur for models having equilibrium states with smaller 
values of the longitudinal to extrinsic propagation velocity,;^ < 1, it is sufficient to 
consider the Newtonian limit 

f 2 <l, c 2 <l. (5.14) 

for which the general eigenvalue equation (4.23) reduces to a form in which the ratio 
X is the only parameter involved appart from the ratio n, the result that is obtained 
being expressible as 

a 3 -4na 2 + (4(n 2 -l) - (n 2 + l)( X 2 - l))a + 2n(n 2 -l)( X 2 - 1) =0. (5.15) 

The results obtained above to the effect that all the roots are real for n = and n = 1 
will still of course apply to this limit in particular, but as soon as we go on to consider 
the next level of excitation as given by n = 2 we find that dynamic instability will 
indeed occur for any state that is "too taut" in the sense of having too small a value 
of x- F° r the equation (5.15) the stipulation that the roots be real is equivalent to 
the requirement that the ratio x be greater than a small but stricty positive critical 
limit value Xc say, i.e. the condition for stability will be satisfied if and only if 

X>Xc (5.16) 

where the critical value is about 1/2, being given approximately by 

Xc^0.47. (5.17) 

Since it is only displacements in the plane of the undisturbed configuration that 
are involved, this limit should be easy to test experimentally by spinning up an 
ordinary elastic string loop such as a thin elastic band on a smooth plane surface. In 
the Hookean low tension limit, x wm be very large so the danger of instability will 
not arise, but as the angular velocity and therefore also the tension of the loop is 
increased one would expect x to decrease monotonically and ultimately to fall below 
the instability threshold (5.17). 

In the relativistic regime relevant to the superconducting string loop problem that 
originally motivated this work the situation is, of course, more complex, depending 
not just on the ratio x but on the absolute values of v and c. A more thorough analytic 
and numerical exploration of the zones of stability and instability in the general case, 
and of their relationship to the predictions of particular superconducting cosmic string 
models, will be postponed for a subsequent article [26]. 
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